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ABSTRACT 



One of the non-linearities which occur in feedback control systems 
is transportation lag. All of the currently popular mathematical tech- 
niques available for the analysis of these systems are severely compli- 
cated by the introduction of transportation lag. The suitability of a 
numerical method employing the Z- transformation to this problem was in- 
vestigated, as well as the particular advantages and disadvantages of 
the other methods. Solutions by the numerical method were compared 
with solutions by other methods and the results obtained presented in 
tabular form. These results demonstrate the superiority of the numeri- 
cal method in determining the transient response of these systems with- 
in reasonable engineering accuracy in a much shorter time. 

The writers wish to express their appreciation to Dr. M. P. Pastel 
for his guidance during the initial stages of the work, and to Dr. G. J. 
Thaler for his advice and suggestions throughout the investigation and 
his aid in the preparation of this paper. 
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1. Introduction 



1. 1 Purpose of the Study 

Feedback control systems are seldom truly linear, but rather 
they are subject to the Introduction of a considerable number of non- 
llnearltles. Some non-llnearltles may be introduced intentionally by 
the system designer. An example of this might be the selection of a 
relay servomechanism with its attendant non-linearity, when a linear 
system could accomplish the task just as well. On the other hand there 
may be no way of avoiding a non-linearity in the system design. An ex- 
ample of this would be the mandatory inclusion of a gear train with its 
attendant backlash. 

A non-linearity which frequently can neither be avoided nor 
neglected is transportation lag. 

Because of the variable influence of transportation lag in 
feedback control systems and its avoidance whenever possible, it has 
received relatively little prominence in the literature. Rather, the 
mathematical tools capable of analyzing systems with transportation lag 
have been indicated without further amplification. While it is true chat 
the tools indicated for the analysis of systems with transportation lag 
will yield satisfactory results, the problem does not end there. These 
standard tools, while yielding accurate results with a reasonable amount 
of work for linear systems, require considerably more work and time when 
the problem is complicated by transportation lag. Consequently, a need 
for a better tool than those currently in use for the solution of pro- 
blems involving transportation lag is indicated. 

In 1954, a paper by Ragazzini and Bergen outlined a numerical 
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methods technique for determining the transient response of linear feed- 
back control systems* The nature of the mathematics which the method 
used Indicated that It might handle the non-linearity of transportation 
lag more conveniently than the aforementioned standard mathematical 
methods generally known to engineers. The purpose of this study Is to 
evaluate the numerical methods technique of Ragazzlnl and Bergen in 
analyzing feedback control systems with transportation lag and to com- 
pare the accuracy attainable and the work Involved In using this method 
with the accuracy attainable and work Involved In using more popular 
methods* 

In doing so. It will be necessary to describe briefly the 
common analytical and graphical methods available for solving linear feed- 
back systems* The Increased difficulty and complexity Introduced by the 
transportation lag Into these methods will be demonstrated* Finally, 
the application of the numerical method to the same problem will be 
presented* 

1*2 Summary of the Results and Conclusions 

The numerical methods technique was compared with the root 
locus, Bode diagram, power series expansion, and second order Fade approxi- 
mation approaches to the problem of system analysis* The Bode diagram 
proved to be the speediest method of all for determining system stability, 
but this method yields no other specific Information regarding the over- 
all system performance* For the full transient response the numerical 
methods technique proved to be faster and more accurate than the second 
order Fade approximation* The root locus method, while retaining Its 
potential for good accuracy, required a prohibitive amount of time and 
tedium to be considered a practical working tool for problems Involving 
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transportation lag. The power series expansion approach to the pro- 
blem proved to be totally impractical for other than extremely simple 
systems. 

In short, the numerical methods technique is the fastest and 
most accurate means of determining the transient response of a feedback 
control system with transportation lag. Moreover, the method can be 
understood and applied by engineers and students whose background includes 
only a first course in linear servomechanisms. In the solution of the 
problem by the numerical methods technique, no more complicated computing 
device than the desk calculator is required. 

In all cases, solutions determined by the various methods were 
compared with the transient response obtained from programming the Fourth 
Order Fade Approximation of the Laplace shift operator on a ten channel 
Donner analog computer. This result was the most correct obtainable 
under the conditions of time and equipment available for the study. 
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Description of the System 



The feedback control system to be studied is shown in block 
form in figure 2-1. In order to carry out the goals of this thesis, a 
number of requirements were set for this system. 

A frequency response diagram crossover of one radian per 
second for all variations of the system was desired* For this reason, 
the value of K in block one will always be equal to one* 

The forward loop transfer function must always be low pass in 
nature* For this reason the following values of G(s) were selected in 
block two* These values also give us convenient first, second, and 
third order systems to investigate* 



G(s) = 1 



2.1 



G(s) 



/ 

O.IS + I 



2.2 



G(s) = 1 

(O.lS-h/)(O.OSs-h0 



2.3 



Block three contains the transportation lag. The nature of 
this block and its analytic representation will be discussed in detail in 
section 4* 



The system will use unity feedback* 
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Figure 2~1, Goneralizec! Block Diagram of Basic System. 
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Linear System Analysis 



3. 1 Discussion 

Various methods of analysis and solution of linear feedback 
systems are available. Two graphical methods currently in use supply 
the necessary information required for the design and analysis of these 
systems. For the purposes of this study, the transient, or time varying 
response characteristics, are more important. 

3.2 Frequency Response 

The Bode network therorems presented a new approach to the 
standard Nyquist stability criterion. One of the great values of the 
Bode diagram is the simplicity with which approximate results of suf- 
ficient accuracy may be obtained. 

The Bode diagram for the linear second order system described 
in section 2 is shown in figure 3-1. Phase margin and gain margin values 
are indicated, and reasonable predictions of system performance may be 
made from these values. This graphical presentation of the open loop 
frequency response is plotted knowing the location of the open loop poles 
and zeros. 

System stability may be inferred from the values of gain margin 
and phase margin on the open loop frequency response. Information from 
the Bode diagram may be transferred to the Nichols chart, and the closed 
loop frequency response may be obtained. Infoirmation from the Nichols 
chart should agree closely with that from a basic Nyquist plot. 

3.3 Root Locus 

A graphical method exists for deteirmining the locus of the 
roots of the closed loop system from the location of the poles and zeros 
of the open loop system. The root locus method essentially demonstrates 
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Figure 3 - 1 . Frequency Response Diagram of Undelayed Second Order System 



the transient response characteristics of the system rather than the 
frequency response. 

The root locus plot for the linear second order system describ- 
ed in section 2 is shown in figure 3-2. Since this method indicates 
directly the effect on the closed loop response of the system gain from 
zero to infinity, it is obvious that this system is stable for all values 
of K. This may also be inferred from inspection of figure 3-1, which 
shows that the phase angle curve always remains above the 180 degree 
phase angle line. 

3*4 Transient Response 

The transient response of a control system contains nearly all 
the parameters required by the engineer for the evaluation of its perfor- 
mance^ A typical transient response is shown in figure 3-3 and the vital 
parameters are indicated. From a study of these values, an estimate of 
the worth of the control system may be made, and the need for compensa- 
tion may be determined. 

The parameters which will be of interest are the followings 
Rise Time - Length of time for the system to first reach a point within 
5% of the desired value. 

Peak Overshoot - Value in percent by which the system overshoots the de- 
sired value. 

Time to Peak Overshoot - Length of time to reach peak overshoot point. 
Settling Time - Length of time for the system to settle within 57. of the 
desired value. 
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4. Transportation Lag 
4. 1 Discussion 



’’Transportation lag” is a terra frequently used in the process 
industry to describe a pure time delay. Other terms used to describe 
this phenomenon include ’’dead time”, ’’delay time”, and ’’true time lag”. 
Figure 4-1 outlines an arrangement which might be used to control the 
thickness of steel plate produced by a rolling mill. A sensor indicates 
the actual plate thickness in voltage form. The difference between its 
output and another voltage representing the desired thickness forms the 
error voltage. The amplified error ultimately positions the rolls. The 
transfer function from plate thickness at the sensor X through the ampli- 
fier, motor, and gear train to roll position, hence thickness at the rolls 
Y, differs in no sense from ordinary transfer functions considered pre- 
viously. From Y to X however a new element is introduced because a finite 
time must elapse before a change in thickness at Y reaches the point of 
measurement at X, d units of length away. If the steel plate proceeds to 
the left with a constant velocity v, any time function of thickness y(t) 
resulting from changes in roll position will be reproduced at the sensor 
without distortion d/v time units later. It is precisely this time quan- 
tum that the term transportation lag describes. 

Transportation lag may be expressed by the following mathe- 
matical relationship: Given f(t) for t ^0, the function f(t-T) differs 

from the function f(t) by the dead time T, where f(t-T) is defined to be 
zero for t^T. 

The setup of problems on computers occasionally requires the 
representation of transportation lag or an expression of the form f(t-T). 
The brute-force approach to the solution of this problem, and one that 
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Figure 4~1* Feedback Control System for Steel Rolling Mill Application 
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has been used in the past, is to generate f(t) in the computer and plot 
f(t) vs. t on a specially equipped recording device. Schemes involving 
such mechanical devices to accomplish the introduction of transportation 
lag have proved to be workable for some problems, but the limitations 
imposed by the physical design of the devices and the introduction of 
human error in their use have made necessary the development of better 
methods. 

The Laplace transform of f(t-T) is 



- sT 

if f(t-T) = 0 for t<T. In equation 4,1 e is commonly called the 
Laplace shift operator. The problem of generating a function f(t-T) 

_ sT 

can therefore be reduced to the problem of generating the function e 
-sT 

Since e cannot be expressed as a rational fraction, the transfer func 

tion for dead time can be simulated only approximately. Thus the remain 

* sT 

ing methods depend on an approximation of the transfer function e by 
an approximate series or technique of network synthesis. 

4. 2 Methods of Simulation 
4.21 Power Series 

X 

The best knovm approximation to the function e is the Taylor 
series expansion 






+ 



iC 

n' 



4.2 



This expansion is convergent for all values of x, but unfortunately 
the rate of convergence is small for x large. For this reason the 
Taylor series is not well suited to the generation of transportation 
lag for problems involving high frequencies or long time constants. 
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This can be seen by replacing s by j u) in the expression for e 

Moreover, in attempting to determine by longhand calculations, the trans 

ient response of other than the simplest system possible involving trans 

portation lag, the introduction of the Taylor series approximation for 
- sT 

e complicates the problem so severely as to render it useless* 

4.22 The Fade Approximation 

In section 4.1 an important theorem of Laplace transforms was 
stated, viz. : 

repeated 4.1 



Thus the transfer function of transportation lag is 






e = 



Lo 

Ei 



Substituting s = j(x) to find the frequency response gives 



= c 






Cos ooT - jrSiNudT 



Thus 



= \jcos^ujT -+ SiN^ixiT - j 



fo 


- \|C05 


Ei 

(t> 


N 

= tAN 



-I 



SW u)T 



00 T 



4,3 



4,4 



4.5 



^ Cos U)T, 

or the frequency transfer function of transportation lag is characteriz- 
ed by an amplitude ratio which is unity for all frequencies and a phase 
angle which varies linearly with frequency over the desired range of in- 
terest. 



A popular method of approximating the transfer function e 
from this point of view is the Fade approximation. It stems from the 
fact that the quotient of any pair of comp lex- conjugate polynomials has 
an amplitude of unity and also possesses considerable phase shift. Ex- 
pressed in its basic form, the Fade approximation is 
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4.6 



X _ UM 

~ (u+^)->0o 



F'u.v (^) 
G-U,V (x) 



where 



Fu,v M 



j/X 

/44 V) /' 



V (v-i) 



-h 



v(v-^) 



(u+v)/u+v-0 Z! 
• z-/ x'" 



+ 



(u+v)(u+ v-i) . . •(U + /) vj 

Gru,v(x) = /— -J±X— ^ U ( U,~/ ) 

(V+u)// (v+u)(v+u.~i) 2 / 



-f 



+ _ U(u-i) ... £./ 

(v-i-u) + (v+\)^! 



4.7 



4.8 



The convergence for this series expansion is quite rapid and often values 

of u and v of 2 will give good accuracy for short time lags and low fre- 

-sT 



quencies. The Fade approximant of e 
tion 4.10 below. 



for u = V = 2 is given in equa- 



-sT ^ i-MsT) + 



4.9 



6 sr 4 /^ 



4.10 



In terms of the original problem, generating f (t-T) when f (t) 
is known, the equation for the second-order Fade approximant can be 
expressed in transfer- function notation as 



f(t-T) ^ D^T^-6or-h /^ 4.1 

f(t) D"T^+4DT4-/a 

The justification of the validity of equation 4.11 is demon 
strated as follows: 









repeated 4.1 
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4.12 



Substituting equation 4.11 in equation 4.1 gives 

J[[-P(t-T)] ^ S^T^- GsT 

llUt)] S^r + (^sT-t/2 

Rewriting equation 4.12 gives 

+^ST+/2)X[f(t-T)J = (s^T^-^sT+ll) 

Performing the inverse Laplace transformation on equation 4.13 gives 

(D'f+^,Dr+l2)f(t-T} = (D'T"-4DT + /2j f(t} 4.14 

Comparison of equations 4.14 and 4.12 reveals that the equations are 
identical. This is, of course, a logical conclusion, as the Laplacian 
operator and the differential operator may be used interchangeably pro- 
viding the initial conditions in a system are identically zero. This is 
true in the system described above. 

-sT 

The fourth order Fade approximant (u = v = 4) for e is 

-sT^ (srf-ZO(sT)'+ leo {ST)^~ sT-^ /4gQ 
^ (5T)U 20{sTf + IS0(5T)\ 8‘iOsTt IGSO 

The Fade approximants may be used for either longhand or com- 
puter solutions. For use in longhand solutions the approximant to be 

-sT 

used is merely substituted for e . It is worthy of note that the 

order of the characteristic equation is raised by the order of the 

-sT 

approximant used to simulate e . If the second order approximant is 
used, two additional roots are Introduced into the characteristic equa- 
tion, and if the fourth order approximant is used, four additional roots 
are introduced into the characteristic equation. 

The technique of programming the Donner Analog Computer for 
handling the transfer function of the fourth order Fade approximant 
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Is included in Appendix Section 9. 

The method of simulating transportation lag by the Fade ap- 
proximants yields what may be termed optimally flat phase characteris- 
tics at the origin but does not fit the phase characteristics as well 
over an extended frequency range. The major drawback to the use of the 
Fade approximants to simulate transportation lag is their response to a 
step input or other input having high-frequency components, since the 
amplitude ratio is always unity but the phase is approximated only for 
lowtuT products. The phase shift characteristics of certain of the Fade 
approximants are shown in Figure 4-2, As a result of the inability of 
these approximants to linearly phase shift high frequencies, the output 
of an approximant ^rings'* badly for high-frequency Inputs, This behavior 
is illustrated in Figure 4-3 which depicts the response of the fourth 
order approximant to inputs of a step and ramp with unity feedback. It 
is noteworthy that the ’’ringing** is considerably less for the ramp in- 
put than for the step input. The reason, of course, is that the ramp in- 
put contains a lower percentage of high-frequency terms than does the step 
input . 

In the solution of problems concerning feedback control systems 
with transportation lag, advantage may be taken of the fact that the typi- 
cal system encountered is, in effect, a low pass filter. Block diagram 
algebra techniques are employed in the solution of these problems. Since 
the commutative law is applicable to block diagram algebra, the actual 
location of the transportation lag transfer function in the forward path 
is immaterial and it may be placed at the end of the path in the computer 
where the high frequency components have been attenuated to the greatest 
possible extent. Frequencies which have been attenuated 40 db, or more 
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UJ T Product (Radians) 




Figure 4-2. Phase Shift Characteristics of the Pade Approximation 
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Figure 4"3(a) . Response of Fovirth Order Fade Approximtion of 
0,9 Second Time Delay to Step Input, 
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Figaure 4"3(b), Response of Fourth Order Fade Approximation of 1.0 

Second Time Delay to Ramp Input TJith Unity Feedback. 
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are considered to be negligible and whether they have been linearly 
phase shifted is of no consequence. This tends to obviate the neces- 
sity of having to linearly phase shift a high uJ T product. Another 
consideration which makes unnecessary the linear phase shift of high uJ T 
products in the solution of the transient response of feedback control 
systems is that high uJ T products mean extremely large phase shifts and 
probable instability. In these cases the transient response is virtually 
meaningless. 

In the great majority of cases the fourth order Fade approxi- 

mant is more than adequate for the simulation of transportation lag. In 

cases where linear phase shift of higher uuT products than those shown in 

figure 4-2 is desirable, higher order Fade approximants or more complex 

- sx 

curve fitting techniques may be used to simulate e . The curve fitting 
techniques are discussed in the literature on analog computers but the 
need for such approaches is questionable. These techniques were devised 
because it was felt the Fade approximants of sixth order and higher fail- 
ed to phase shift low uuT products accurately. In the course of this 
study the phase shift characteristics of the sixth and eighth order Fade 
approximants were calculated. The results are shown in figure 4-4, On 
the basis of these results it may be inferred that the order of the Fade 
approximation may be increased to increase the magnitude of the ojt pro- 
duct that is linearly phase shifted without distorting the phase shift of 
low (jUT products. Because the Fade approximation may be expressed in a 
general form from which any order approximant may be developed, it is 
felt that the use of curve fitting techniques requires a needless expendi- 
ture of effort without any tangible benefit. 
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Figure 4-4. Phase Shift Characteristics 
of the Fade Approximation 
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5. 



Effect of Time Delay on Linear System Analysis 



The addition of the time delay into the linear control system 
under investigation introduces several factors which tend to make the 
analysis of the system more difficult. Because the time delay may basical- 
ly be represented by an infinite series of terms, this creates an infinite 
number of roots. The effect of these complications on the standard graph- 
ical analysts techniques will now be investigated. In all cases, the 
system with a time delay of 1.0 seconds will be demonstrated. 

5. 1 The Frequency Response 

The standard Bode diagram is shown in figure 5-1 when the time 
delay is included in the system. The difference between this plot and 
figure 3-1 is readily apparent. The phase margin has been reduced by a 
considerable amount. The gain margin is approximately 3.3 db which 
compares with the undelayed system which was stable at all gains. 

The time delay has introduced additional branches into the dia- 
gram. These branches appear due to the periodic nature of the time de- 
lay, and the increasing phase shift introduced at higher values of fre- 
quency. It is evident that the gain margin for each of the branches in- 
creases. This is analogous to the statement that the roots of the pri- 
mary branch would go unstable first if the gain were raised. 

The construction of this plot is simple and straightforward. 

For the system in question the open loop transfer function is 



The gain vs log uj curve is identical to the undelayed system 
curve since the time delay introduces no additional gain into the system. 
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Figure 5“1« rVeqxiency Response Diagram of Second Order System 
\7ith Transportation Lag of One Second. 



The phase angle vs logo) curve is drastically changed since 
the equation for the phase angle is 

(j) = ^0 -f ARC 0./ Cj -f CU 7 d 5^2 

It is the presence of the additional term of the frequency times 
the time delay which introduces the additional phase shift and the extra 
branches. 

5.2 Root Locus 

The identical second order system with a transportation lag 
of 1.0 seconds also appears greatly changed when investigated by a 
root locus plot. 

Figure 3-2 showed the undelayed second order system which is 
stable for all gains. Figure 5-2 has the additional construction de- 
tail which we will require to obtain the root locus. In the same manner 
as the minus 180 degree phase shift line was obtained, lines of constant 
phase shift of 170, 160, 150, 140 degrees, and so forth, may be plotted. 

In addition, lines of constant phase shift caused by the one second time 
delay are also shown. In the non-dimensionalized plot shown, the time de- 
lay will produce a complete 360 degree phase shift every 0.628 radians. 

The resulting root locus is shown in figure 5-3. It is plot- 
ted by combining points of the basic system and the time delay which 
total to a 180 degree phase shift. It is evident that this system will 
go unstable with only a small increase in gain, and as a matter of fact, 
the gain at the stability limit just before the locus crosses the axis is 
approximately 1.48, the same value of gain margin obtained from the Bode 
diagram. 
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Figure 5-3. Root Locus of Second Order System T(ith 
Transportation Lag of One Second. 
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An extremely interesting result of this plot is that upper 
branches of the root locus appear. Only two of these are shown although 
an infinite number exist. Here again they are a result of the 360 degree 
periodicity introduced by the time lag. It is also noted that at the 
same gain, the additional roots on the upper branches are attenuated more 
than the root on the primary branch. This substantiates the Bode diagram 
findings, which show that the roots on the primary branch would go unstable 
first. 



A continuation of the graphical solution described above for 
other values of time delay will yield other root loci curves similar to 
the one in figure 5-3. 

The gain at the stability limit for each time delay was calculat- 
ed. The product of the time constant of the system and the gain, K , 

was plotted against the time delay on logarithmic cross-section paper. 

The result is a straight line as shown in figure 5-4. It was therefore 
possible to write a single equation which will describe the gain at the 
stability limit in terms of the time constant of the second order system 
and the time delay. 






0 . 

.s? 



5.3 



r A To 

Thus for second order systems if two of the above parameters are 
known, the maximum value that the third parameter may have can readily be 
calculated. 

5.3 Transient Response 

Neither of the graphical methods discussed here will directly 
produce the transient response of any system. The Bode diagram will only 



27 



I«IM' II i — t fi»l^ * 

*’^"'mi'maillh-^lmmm Ahm> «■ » *J W«b» • *^' 

■■■^ i4ii0l 
If • m #MA|jHitlM^* I 




II £««» fMlCMV «l iKi II IflillllMil ^1 «MM||| 





Figure 5-4. Non-Dimens ionalized Gain at Stability Limit vs. 

Transportation Lag for Second Order System 
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yield values from which an approximate transient response may be infer- 
red. The root locus method provides the roots of the system on the s 
plane, but then the results must be transformed to the time domain and the 
equation of the output determined. The transient response must be obtain- 
ed by substituting values into this equation and obtaining the output at 
each specific time. Complicating this procedure is the fact that an in- 
finite number of roots exist on the root locus, and enough of these roots 
must be considered to obtain reasonable accuracy. Each time an addition- 
al set of these complex conjugate roots is considered, additional time and 
complexity are added to the above procedure. 

The transient response of this delayed second order system is 
shown in figure 5-5. It is here that the influence of the time delay on 
the second order system is most apparent. The original second order 
system was an overdamped system. It had no overshoot. The time delay has 
deteriorated the system to the point where the overshoot is approximately 
70 percent. 

It is interesting to note that a time delay of about 1.4 seconds 
would cause this system to be completely unstable. 

The effect of the time delay on the transient response there- 
fore is to increase overshoot and settling time. It can make a normally 
overdamped system highly oscillatory with a delay of only one second in- 
troduced. 
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Transient Response by Numerical Methods 



6. 1 Discussion 

A mathematical technique now exists for the direct solution of 
the output of a linear continuous feedback system. An extension of this 
numerical method to the solution of systems with transportation lag is 
easily accomplished. 

6.2 The Z - Transformation 

A continuous time function f(t) is sampled by means of a sampl- 
ing switch S at equal sampling intervals of time T in figure 6-1. The 
output of the sampler f^(t) is a sequence of pulses which may be represent- 
ed for mathematical purposes as a series of impulses or delta functions 
whose areas are equal to the amplitude of f(t) at the respective sampl- 
ing instants. Thus 

rw= ^(t-hT) 6.1 

B — cO 

The Laplace transform of the sequence of impulses is given as 

l4(nT) 

The Laplace transform of such pulsed functions, therefore, is in 
— nTs 

terms of e , Replacing this exponential term by a new variable Z, where 
= e accounts for the name Z- transform. One useful characteristic 
of the Z-transform is that the infinite series in equation 6.2 may usual- 
ly be expressed in closed form so that tables of Z-transforms may be de- 
veloped in the same manner as tables of ordinary Laplace transforms. 

One of the big advantages of the Z-transform is that it is possible 
to obtain the pulsed output of a linear system by expressing a pulse trans- 
fer function which relates the input and output Z- transforms. The pulse 
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f*(t) 




Figure 6-1. Sampled Time Function 
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transfer functions can be expressed in closed form and tables of Z-trans- 
forras are available in most standard texts covering sampled-data systems. 

An abbreviated table is given in Table 6-1. 

6*3 Application to Linear Continuous Systems 

The application of the numerical method using Z-transform 
technique will be illustrated in a linear continuous feedback system as 
shown in figure 6-2 (a). The system is mathematically converted into a 
sampled system by placing synchronized samplers at input and output and 
at point X in the feedback path in figure 6-2 (b). The sampling frequency 
should be high enough to cause the error introduced by the sampled system 
to be small. By locating the sampler in the feedback loop and requiring 
that the transfer function G(s) be low pass in nature, then only the for- 
ward gain of the original system need be considered. Later examples will 
show the effect of the ratio of sampling frequency to the Bode cross- 
over frequency u)^ of the original system. 

Following the sampler S in the feedback loop is the element B(s), 
whose function is to reconstruct the original output of the forward path 
as faithfully as possible. This element reconstructs a polygonal approxi- 
mation to the actual function by converting each pulse into a triangular 
generating function as shown in figure 6-3. This requires the response of 
B to a short pulse or impulse of the same area to be a triangle as shown 
in figure 6-4(a). This response is not physically attainable since it 
requires a slope at negative time, but this has no effect on the purely 
mathematical solution. 

The combination of sampling switch operating with a period T 
and the triangle function generator B replace the direct connection 
which existed in the original system. The errors introduced Into the 
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Table 6-/ 



Time Function Laplace 7ran>?form Rs) Z- IransfoRki 



S f t ~ hT) 
U.(t) 
t 

ii! 

2 

.-at 



l-C 



-at 



-nTs 



5 

/ 

f 



S+Ck 

q 

S(s+c\) 



-h 



l-Z 



-/ 






,-/ 



-/\2. 



(i-rn 

2(1-2'')^ 



-aT^-/ 






z~' 



S//V qt 



a 






2~''S‘/n gT 



I -22 'cos aJ + 2' 



Cos qt 






/~2~cos gT 
-£Z cos 0 .T + E 
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(b) 



Figure 6-2. (a) Continuous Feedback System 

(b) Sampled Model for Computation 
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Figure 6-3. Polygonal Approximation of a Function f(t) 





Figure 6-4. (a) Generating Triangle 

(b) Ramps I'/hich Produce Triangle 
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solution are a result of the Inability of this combination to faithfully 
reproduce the feedback signal. 



The transfer function B(s) of the triangle function generator is 
obtained from the Laplace transform of the triangle function in figure 
6-4(a). This triangle is a result of the combination of three ramp func- 
tions as shown in figure 6-4(b). The Laplace transform of these three 
functions is; ' TS p I _-T5 



rs" TS" 



T5 



z e 



or 









6.3 



The relation between the Z-transform of the output and input in 
figure 6-2 (b) is given by 

0/(2) zz 6.4 

I + 

and yields the pulse transfer functions for which tables are available. 

If the response of a linear system to an input F(s) is desired, 
first replace it with the system shown in figure 6-2 (b). A proper choice 
of sampling frequency (jLj^ is made. The pulse transfer functions RG*(z) 
and HBG*(z) are computed, and combined according to the basic laws of 
algebra, yielding a form of a ratio of polynomials in powers of z. To 
obtain the output pulse sequence in the time domain, a simple process of 
long division is required. The division process will result in a power 
series of the form 



The magnitude of the output in the time domain is just the coefficient 
of each of the terms in this series at a time corresponding to the 
power of z. Thus the output at the fourth sampling instant is a^ in the 



expansion. 
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The long division process is continued until the desired 



number of points have been obtained. The effect of the sampling period 

T , or of the sampling frequency will be elaborated on in section 6.5. 
s 

6.4 Application to Linear Continuous Systems with Transportation Lag . 

It has been shown in section 4.1 that the transfer function of a 

-Ts 

transportation lag or simple dead time lag is merely e where T is the 
value T Q , or the time lag. 

If this transfer function is inserted into the linear system of 
figure 6-2 (b) as Gj(s) then the solution for the output will become 




1 + 



6.6 



The value of this numerical method employing pulse transforms 

-T s 

is now easily seen. The z-transform of the time delay e ^ simply be- 
comes 2 ^ from the basic definition of the z-transform itself. Since the 
remaining terms of fl*(z) are also polynomials in powers of z, this 
does not add any more complexity to the basic long division solution of the 
transient response. 

It must be noted that the T introduced in the polynomial recon- 
struction and the z-trans forms in Table 6-1 is T^, the sampling interval. 
This must be distinguished and kept separate from T q ,the time delay. 

One of the basic rules for the method is that the sampling time 
T^ must be equal to or an integral sub-multiple of the time delay T^ . 

It will be necessary to express 6^(z) only in terms of the sampling time 
T^ in order that the value of the transient response will be obtained at 
each sampling interval. This method will be shown in the sample solu- 
tions in section 7. 
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6. 5 Requirements for Reasonable Accuracy 



The ratio of to the crossover frequency u-?q of the open 

loop transfer function will determine the relative accuracy of the system. 

The comparison of sampling period T with the expected settling time and/ 

s 

or time delay in the system will determine whether the method will supply 
the required answer within a reasonable number of long division steps. 

Ideally, to obtain extremely good accuracy (within 1% of known 
answer) it is necessary to select a value of u>^ which is twice the fre- 
quency at which the open loop system is down 40 db on the Bode diagram. 

For the systems investigated here and explained in section 2, is approxi- 
mately 1 rad/sec. For the system where G(s) = ~ , this would require 

s 

=s 200 rad/sec. Since T^ = ^ this would require a T^ of approximately 

0,03 sec. With a normal settling time for this system of 4 seconds this 

4 

would require a total of — or 133 long division steps to obtain the 

, v/ ^ 

transient value at the end of 4 seconds. In order to reduce the computa- 
tional work and to remain within approximately an accuracy of 5%, vari- 
ous ratios of (J, to were selected and results tabulated against re- 
suits of former traditional methods of solution. 

The problem of reducing computational work is further aggravated 
by the introduction of the time delay into the system. One basic re- 
quirement of this numerical method is that the sampling time T^ must be 
equal to or an integral sub-multiple of the time delay T ^ , For systems 

with a long settling time, for example 5 seconds, and a short time delay 
of 0,10 seconds, it would require 50 long division steps to produce the 
value of the transient response at 5 seconds. This is only if it would 
be possible to select T^ = T ^ , In the event that the value of T^ must 

be chosen smaller to satisfy the requirement of the ratio to ui^ for 
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accuracy, then it would require 50n long division steps, where 

n ^ 

It will be demonstrated in the next section that this method 

will provide reasonable engineering accuracy (approximately 5%) with a 

small investment of time. This can be accomplished with a ratio o f LO^ 

to as low as possible, and a T which Is the lowest sub-multiple of 

s 

T ^ as long as it is not greater than . 

This method can give extremely good accuracy at the expense of 
an extremely long and repetitious series of calculations. Its advantage 
lies in the fact that where the ratios mentioned above are favorable, it 
should be possible to compute the transient response of a system in a 
short period of time with no more powerful tools than algebra and a slide 
rule or desk calculator. 

The next section will demonstrate the accuracy and simplicity 
of the method in solving simple feedback systems. Computational time 
and labor will be compared with other methods currently available to the 
engineer. 
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7 , Comparison of Numerical Method with Other Methods of Solution 
?• 1 Numerical Method Solution 

The numerical method solution to the second order system des- 
cribed in section 2 will now be demonstrated in detail. The system to 
be solved is shown in figure 6-2 (b). The input is a unit step. The time 
delay will be one second. The sampling time interval will also be one 
second. The table of results in section 7.5 will also show the solution 
results for other sampling intervals. 

The output in the Laplace notation is 

-S-TL 



where 



= 






/o 



s^(s+io:> 



/ + 



S+/0) 



7.1 



IIS^ 



/ - sH \ 2 

(/- e ) 



7.2 



Combining 



0c(^) = 



IOC 






(S-tlO) 



' (S+/o) 



7.3 



It is important to note that the T in the time delay represents 

Td , the delay time, whereas the T in the expression for B(s) and also in 

the Z transforms represents T , the sampling time. In the algebraic work 

s 

in combining the expressions it will be necessary to express Tt> in terms 

of T and obtain the output equation in terms of T only. In this ex- 
s s 

ample, since To equals T equals one second, it will not be necessary to 

s 

account for this difference. 

By expressing numerator and denominator separately in partial 
fraction expansions, and then transforming each separately by means of 
the transforms in table I of section 6 we obtain 
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7.4 



lo(l~^"Y 

Zoa-//s>2-'^ //C,2'2- + 22'^ 

20O( /-2-'J) 

Combining numerator and denominator results in 

0 .(£) - l 20 Z -^-^ lor ^_ 7.5 

200 - J/^> + 23^f 2'2_ // Y g-'' 

By carrying out the long division indicated in equation 7.5 
we obtain a power series in Z 

6c iSZS-Z-"^ ^0.9/ 2."^+ • • • 

We may directly transform this expression into the time domain. 
It is not really necessary since the only effect on the output will be 
to express the series in terras of the sampling periods. In either case 
the coefficients are the same and are the value of the transient response 
at each sampling instant. 

^ I.S3l(t-3T)-t J.3SS(t’HT)i- • • - 7.7 

The introduction of a different time delay or sampling time 

into the problem merely changes the algebraic relationship between 

and T and once this accounting has been accomplished, the algebra and 
8 

final solution work is similar to that shown here. No additional diffi- 
culty is introduced into the problem. 

7. 2 Power Series 
7.21 First Order System 

The solution of the transient response of a first order system 
with transportation lag by the employment of a power series expansion is 
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quite straight forward, although tedious and time consuming. The 
system under consideration is that shown in figure 7- 1(a) where T is 
equal to one second. 



The closed loop transfer function may be written 

&c _ S 

0. - 

For a unit step input 






s 



Then 



0r ^ 






/ + 



•Ts 






rT5\-l 



0 - 



-Ts - 1 

The quantity ( 1 + e ) may be expanded in a power 



with the following result. 



6( 



_ e 



-Ts 



■ g-.Ts g-3Ts 

L' s S'- S' 






0c = 






5 



^-3Ts 

7rn~ r-v '<r 






Since T = 1 



-s 



a- 






C- 3 I OH 












+ 



Taking the inverse Laplace transform of each term yields 

0cCt) " ^ +A(t-3j^ 



7.8 



7.9 



7,10 

series 



7.11 



7.12 



7.13 



7.14 
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Figure 7-l(a) , First Order Feedback Control System 
V.'ith Transportation Lag. 




Figure 7-l(b) , Second Order Feedback Control System 
VJith Transportation Lag. 
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The transient response may now be determined by substituting 
values of t into equation 7.14 . The magnitude of the output at any 

time t may be determined as exactly as desired by evaluating all terms 
and by evaluating each term to a sufficient number of decimal places. 

The amount of time required to obtain a transient solution by this method 
increases directly as the time interval between sampling instants is de- 
creased. It also increases as the settling time increases since increased 
settling time increases the number of terms which must be evaluated at 
higher values of t» The sample problems solved by this method required 
about twelve hours to complete. 

7.22 Second Order System 

While the use of a power series expansion is a straight forward 
method of determining the transient response of a first order system with 
transportation lag, the method requires far too much work to be considered 
practical for higher order systems. The following example will demons- 
trate why the method fails ^as a practical working tool. The system under 
consideration is that shown in figure 7-l(b) where T is equal to one 
second. 



The closed loop transfer function may be written 



Then 



-5T 






For a unit step input 
0-R. - 



^(O./S + l ) 



I i- 



S 



,-sT 



9c- 



.-ST- 






S"^(o./s*i) 



,-sT 



/ + 



=-ST 



S(o-^S4i) 



-sT \-l 






V S(ais+i) 



7.15 



7,16 
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The quantity 



may be expanded in a 



/ + 



>-rT 



power series with the following results. 

0c = e’"'' ^ 



S*(0.(S+4^ 

^-3ST 



S(pj 

g results 

r ^-sT ^-2sr 

I / — — ^ -f — ^ 

[_ S(d.lS+i) sXows-41)' 



-1 

.'/S-Hr J 



+ 



7.17 



s*(a/s+()^ 

At this point it may readily be seen that the task of determin- 
ing the inverse Laplace transform of each term is momentus and the work 
involved increases tremendously with each increase in the order of the 
term. Thus it is maintained that the use of a power series is impracti- 
cal in determining the transient response of systems of second order or 
higher which are complicated by transportation lag. 

7. 3 Fade Approximation (Arithmetic ) 

The transient response will now be determined for the system 
of section 7.1 by simulating the Laplace shift operator with the second 
order Fade approximant. The major steps in the solution will be outlined 
but the tedious arithmetical calculations will of necessity not be repeated 
in their entirety. 

The system in block diagram form is shown in Figure 7- 1(b). 

Since the duration of the transportation lag is one second the Laplace 

- s 

shift operator becomes e . As indicated in section 4.22 






5 - fcS -f IZ 

+ IZ 



7,18 



Substituting the Fade approximant in for the Laplace shift 
operator and manipulating by block diagram algebra yields the closed 
loop transfer function for the system, 

0 



0R 



/ + IZ') 

I S (5+(oXsV4S+(2.) 



7.19 
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For a step input, 



/ 



7.20 



n - _i / lOS^-^OS -I- 12.0 

5 \ Eis^+Cos + /2o 

The transient response of this system to a step input is con- 
tained in Laplace notation above. The problem of reverting to the time 
domain remains. 

The first step is the determination of the four roots of 
4 3 2 

S + 16S + 82s + 60S + 120. Since these roots are two complex conju- 

gate pairs whose real and imaginary parts are not necessarily integral 

the Lin method was selected for their solution. The first trial divisor is 
2 2 

82S + 60S -f 120 which equals S + .733S 1.46. Five iterations were 

2 2 

necessary to yield factors of ( S + .463S + 1.65 ) ( S + 15.548^ 73.2 ). 
Application of the quadratic formula further reduced the factors to ( S + 
.232 + jl.264 ) ( S + 7.77 + J3.59 ). 

Next, the method of residues was employed to determine the co- 
efficient of the inverse transform of each of the five poles. 

Thus 






9 , (t)= I* (-. 3/2 (-.mijMs) 



(-7.V±i3.S^)t 

t ^ 



7.21 



This further reduces to 

9c(t)= {z cosl- 2 (-Ht^ 3 .H(^ SIN !.2(^H±) 

cos 3S^t t 2.63 siNB.S^t) 

The full transient response of the system may now be plotted by 
evaluating 0c (*^) various values of t. 
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The time required for the solution of the transient response 
of the system by this method is between seven and eight hours. It is 
worthy of note that the roots of the characteristic equation determined 
by this method are reasonably close to the values obtained from the more 
accurate but more lengthy root locus solution. Arithmetic solutions 
making use of the Fade approximation were limited to the use of the 
second and third order approximants because the time required to solve 
a problem when higher order approximants are used becomes unreasonably 
long. 

7.4 Fade Approximation (Computer ) 

In section 4.22 it was pointed out that the fourth order Fade 
approximant was more than adequate for handling the vast majority of prob- 
lems involving transportation lag. Accordingly, all problems solved by 

arithmetical techniques were also solved on the Donner Analog Computer 

-sT 

using the fourth order Fade approximant to simulate e . The solution 
to the problem by this method was established as the control against 
which the arithmetic solutions were compared. The fourth order approxi- 
mant was chosen because it occupied only six amplifiers on the Donner 
Computer thus leaving four amplifiers free on which the remainder of the 
system could be programmed. The computer program for the fourth order 
Fade approximant is contained in the Appendix. This program was specifi- 
cally written for a range of values of transportation lag between 0.1 
seconds to 1.0 seconds. In each amplifier only one resistor had to 

be changed to change the value of the transportation lag between the limits 
stated above. For transportation lags outside of these limits the program 
should be modified so that does not exceed the bounds of 0.1 to 1.0 
megohm. 
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7 • 5 Table of Results 



This table of results Is an intent to show the relative ac- 
curacy of the various methods, and to compare the solution time in- 
volved to achieve this accuracy. The first order system described in 
section 2 was used to establish the basic comparisons for accuracy, and 
to demonstrate the feasibility of the numerical method. Results indi- 
cate that the numerical method approach provides results of comparable 
accuracy to other methods with substantial savings in computational time. 
Results also indicate that increased accuracy may be obtained with the 
numerical methods solution by increasing the sampling frequency. It 
must be pointed out, however, that the time to arrive at the solution 
increases directly in proportion to the increase in the sampling fre- 
quency. 

Table 7-1 contains the basic information from which general 
conclusions about relative accuracy will be based. The results indicate 
the value of the transient response for a first order system with a one 
second time delay as obtained by various methods. Results are shown at 
one second intervals. The power series solution is the most accurate 
available and will be taken as the standard by which the other methods 
will be compared. 

Table 7-2 contains similar information for the first order sys- 
tem with a 0,5 second delay. The numerical method solution for a ratio 
of GJs to ratio of 12.56 is compared because it achieves the stated 
requirements of the study of less than 5% deviation. It will be the ratio 
used in higher order systems. 

Tables 7-3 and 7-4 display results for the second and third order 
systems respectively for a 0.5 second delay. The basis for comparison in 
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these systems was the computer solution employing the fourth order Fade 
approximant. The deviations here from the standard are of the same 
orders of magnitude as experienced in the first order system. In these 
results, a plus deviation indicates a result higher than the standard, 
and a minus deviation a result lower. The numerical method solution 
when compared with the standard characteristically maintains the same 
sign at the same relative point of the transient response. That is, the 
transient response as determined by the numerical method coincides with 
or leads the standard, and would therefore tend to show a slightly short- 
er rise time and time to peak overshoot. Its accuracy in determining the 
peak overshoot and settling time was well within the engineering accuracy 
prescribed for this investigation. 

It must be understood that the solutions by the numerical methods 
technique yield ’’exact** values only at the sampling instants. Interpola- 
tion may be done to determine information at times between these sampling 
instants with assurance of reasonable accuracy as the magnitude and alter- 
nating signs of the deviations of the tabulated results from the standard 
indicate. 

All solutions, both computer and analytic, employing the Fade 

approximation exhibit large deviations at the first few tabulated times. 

These are due to the ’’ringing” discussed in section 4,22 and illustrated 

in figure 4-3, This is an unavoidable disadvantage of the Fade approxi- 
- sT 

mat ion of e 

It cannot be over emphasized that the time required to obtain 
the transient solution by the numerical method was in most cases from 
one-fifth to one-tenth the time required by previously used analytical 
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and graphical methods. Its value also lies in the fact that the pro- 
cedure is simple and easy to follow. Once the engineer or student 
becomes familiar with the procedures, it is possible to recognize and 
correct computational errors before they are compounded in following 
steps. This is simply not possible in the other analytical methods 
until the final step of the solution. Similarly, a change in the 
transportation lag or the sampling interval is easily handled, and 
a great deal of the algebraic manipulations need be done only once 
for each basic system. 
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TABLE 7-1 



COMPARISON OF TRANSIENT RESPONSE AS OBTAINED BY 
POWER SERIES WITH OTHER ARITOMETIC METHODS 
FIRST ORDER SYSTEM - TRANSPORTATION LAG EQUAL TO ONE SECOND 



ime 

Sec, 


Power 

Series 


2nd 

Order 

CArlth.) 


7 , 

Dev. 


Num. 

Meth. 

^'-6.28 


7 . 

Dev 


Num. 

,, Meth. 
Sm2.56 


7 , 

Dev. , . 

_<uOs 

aJo‘ 


Num. 

Meth, 

=25.12 


7 , 

Dev. 


0 


0 


.011 


- 


0 


- 


0 


- 


0 


- 


1 


0 


.291 


- 


0 




0 


- 


0 


- 


2 


1.00 


1.013 


+1.3 


1.00 


0 


1.00 


0 


1.00 


0 


3 


1.50 


1.485 


-1.0 


1.50 


0 


1.50 


0 


1.50 


0 


4 


1.167 


1.158 


-0.8 


1.250 


+7.1 


1.1875 


+1.8 


1.1718 


+0.41 


5 


.791 


.795 


+0.5 


.875 


+10.6 


.8125 


+2.7 


.7969 


+0.75 


6 


.8416 


.848 


+0.8 


.8125 


-3.5 


.8320 


-1.1 


.8392 


-0.28 


7 


1.057 


1.055 


-0.2 


.9688 


-8.3 


1.0331 


-2.3 


1.0509 


-0.58 


8 


1.100 


1.098 


-0.2 


1.0782 


-2.0 


1.0993 


-0.1 


1.1022 


+0.2 


9 


1.004 


1.005 


+0.1 


1.0548 


+4.9 


1.0205 


+1.6 


1.0087 


+0,37 


10 


.948 


.950 


+0.2 


.9884 


+4.2 


.9564 


+0.9 


.9543 


+0.66 


11 


.9798 


.981 


+0.1 


.9669 


-1.3 


.9719 


-0.8 


.9916 


+1.20 


12 


1.021 


1.019 


-0.2 


.9894 


-3.1 


1.0108 


-1.0 


1.0329 


+1.15 
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TABLE 7-2 



COMPARISON OF TRANSIENT RESPONSE AS OBTAINED 
BY POWER SERIES WITH NUMERICAL METHOD RESULT 
FIRST ORDER SYSTEM-TRANSPORTATION LAG EQUAL TO 0.5 SECOND 



Time 

Sec 


Power 

Series 


Numerical 

Method 


% 

Dev. 


0 


0 


0 


0 


0.5 


0 


0 


0 


1.0 


0.5 


0.5 


0 


1.5 


.875 


.875 


0 


2.0 


1.028 


1.031 


+0.29 


2.5 


1.039 


1.055 


+1.06 


3.0 


1.021 


1.033 


+1.18 


3.5 


1.006 


1.011 


+0.5 


4.0 


.999 


1.000 


+0.1 
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TABLE 7-3 



COMPARISON OF TRANSIENT RESPONSE FROM COMPUTER SOLUTION 
WITH NUMERICAL METHOD AND PADE APPROXIMANT (ARITHMETIC) 
SECOND ORDER SYSTEM - TRANSPORTATION LAG EQUAL TO 0.5 SECOND 



Time 

Sec 


Computer Numerical 

Method 

12.56 


% 

Dev 


Pade 

Approx. 

(Arith) 


7, 

Dev. 


0 


0 


0 


- 


.01 


- 


0.5 


0 


0 


- 


.0018 


- 


1.0 


.31 


.35 


+12.8 


.404 


+30.4 


1.5 


.77 


.803 


4 4.3 


.846 


+ 9.9 


2.0 


1.04 


1.078 


+ 3.7 


1.073 


+ 3.2 


2.5 


1.115 


1.153 


+ 3.4 


1.111 


- 0.5 


3.0 


1.100 


1.110 


+ 0.9 


1.066 


- 3.1 


3.5 


1.050 


1.037 


- 1.2 


1.017 


- 3.1 


4.0 


1.015 


.984 


- 3.1 


.992 


- 2.3 
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TABLE 7-4 



COMPARISON OF TRANSIENT RESPONSE FROM COMPUTER SOLUTION 
WITH NUMERICAL METHOD AND PADE APPROXIMANT (ARITHMETIC) 
THIRD ORDER SYSTEM - TRANSPORTATION LAG EQUAL TO 0.5 SECOND 



Time 

Sec. 


Numerical 

Computer Method 

^^=12.56 


% 

Dev. 


Pade 

Approx. 

(Arith) 


7, 

Dev. 


0 


0 


0 


- 


.146 


- 


0.5 


0 


0 


- 


.0203 


- 


1.0 


.26 


.35 


+34.6 


.359 


+38.0 


1.5 


.73 


.803 


+ 10.0 


.821 


+12.5 


2.0 


1.050 


1.078 


+ 2.7 


1.09 


+ 3.8 


2.5 


1.150 


1.153 


+ 0.3 


1.150 


0 


3.0 


1.115 


1.110 


- 0.5 


1.099 


- 1.4 


3.5 


1.060 


1.037 


- 2.2 


1.029 


- 2.9 


4.0 


1.005 


.984 


- 2.1 


.987 


- 1.8 
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8 • Conclusions and Recommendations for Further Investigation 



The investigation into the effect of transportation lag on 
feedback control systems reported on in this thesis has demonstrated 
four ideas of practical significance. It has also opened up a number 
of avenues for further investigation. 

The primary conclusion which the investigation supported is 
the applicability of the numerical method technique devised by RAGAZZINI 
and BERGEN to the problem of determining the transient response of linear 
feedback control systems which are complicated by transportation lag. Of 
all the mathematical and graphical tools available to the engineer it 
offers the optimum combination of speed, accuracy, and simplicity. It 
may be used by students whose background includes courses only as advanc- 
ed as linear feedback control systems, provided the Z-transform is included. 

The deviation of the transient response as determined by the 
numerical method from that determined by the techniques currently in use 
may be held to less than three percent. Meanwhile the time required to 
determine the transient response has been reduced by more than seventy- 
five percent. Thus the decision of whether compensation is required or 
not may be made much more rapidly than before. The form the compensation 
will take will, of course, vary from problem to problem and from designer 
to designer. 

The suitability of the numerical method to the particular pro- 
blem of transportation lag suggests an area of further investigation. In 
determining the transient response the system has been mathematically con- 
verted to a sampled-data system, and is described by Z-transforms. In 
order to deteirmine whether the system will meet specified performance 
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criteria it is necessary to calculate the transient response*, If the 
system is not satisfactory, compensation techniques applicable to con- 
tinuous feedback systems would be employed and the transient response of 
the compensated system recalculated. However, since the system is al- 
ready a mathematical sampled-data system, it would appear that the suit- 
ability of sampled-data system performance criteria and compensation 
techniques to the problem is worthy of investigation. Should these 
techniques prove successful the reverting back and forth from the Laplace 
transform to the Z-transform for compensation could be avoided and a 
greater saving in time and tedium effected. 

The second noteworthy conclusion of this Investigation is 
that the Bode diagram not only offers the fastest tool for the deter- 
mination of system stability but it also will indicate which set of roots 
will go unstable first. This information is made available by no other 
technique except the root locus and a conservative estimate would place 
the time differential using the root locus at a full order of magnitude 
higher than the Bode. 

The criterion of dominance of a set of roots in a higher order 
system provides the basis for the extension of the performance character- 
istics of a second order system to these higher order systems. From the 
information in the Bode and root locus graphical plots, it may be practical 
to construct a set of second order curves with transportation lag which 
might have applicability to the analysis of higher order systemSp 

A third significant conclusion is that the Fade approximation 
for the Laplace shift operator is considerably more linear in its phase 
shift characteristics than the recent analog computer literature would 
indicate. An inspection of figure 4-4 will reveal that the Fade approxi- 
mation provides a completely linear phase shift until that point where It 
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deviates from the ideal. Thus, the simulation of time delay on the analog 
computer requires merely the choice of a high enough order Fade approxi- 
mant which will linearly phase shift the frequencies desired. It would 
appear that the curve fitting techniques currently in existence are not 
justified due to their increased complexity and lack of a simple repeti- 
tive scheme for programming, 

A fourth interesting result of this study is the information on 
stability limits for the second order system, I^en plotted on logarithmic 
cross-section paper as shown in Figure 5-4, the non-dimensionalized gain 
vs. the time delay produced a straight line curve from which it was possible 
to write an equation describing the relationship between gain at stability 
limit and transportation lag for any second order system. 

This result is obtainable from root locus or frequency response plots 
and is easily presented because of the ease in non-dimensionalizlng a 
second order system. Here again, a further extension of this informa- 
tion would be possible for higher order systems. The method of display- 
ing the information would depend on the way in which the additional time 
constants of the higher order systems were handled. 
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9 , Appendix - Computer Program for Fourth Order Fade Approximant 

The Fade approximation of is discussed in Chapter 4, An 
analog computer program for the fourth order approximant which was used 
in the study is developed here. 

The Fade approximation states that 



X _ L/M 



G'u.v (^) 



9.1 



where 






VX 



(u+v) I L 



+ 



(V-l) 



(u-vv)(u-^v-i) 2 ! 



, + v(v-\)--a-l X'' 

(u+v)(u+v-l) • •- (u+i) vj 



9.2 



and 



&u,v(^) = / - 



ux 



u/u-i) 






(-<r 



u.(u-0 



• Z-l X 



(V+U^(V+Lt.-l) . . *(V + l) Ul 
Substituting u = v = 4, and x = -sT yields 

-5T _ (STf - 20(sT)^-f g^O(sT) + i(ogO 

^ ' (STJ%20(sT)"4/FO(sT)-4?^0(5T} + /G?O 



9.3 



9.4 



e 



-sT 



20 . /XO _ SHO 

I - ST {ST)^ (STY 



4 



I (oSO 
(ST)^ 



h 



Zo 



/So 



ST (STf 



S40 

(ST)= 



-4 



/6SO 



As shown in Section 4.22 



-sT 






9.5 



9.6 
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Since the differential operator and Laplace operator are inter- 
changeable when initial conditions are zero, equations 9.5 and 9.6 may be 
combined as follows 

/ -h - g'Z-O , I (o 20 

! DT 

I 4- ^ _1_^0 /480 

' DT D*-T^ D‘*T*' 




The final form of equation 9.7 from which the analog computer program 
shown in figure 9-1 was developed is 

((t-T) = kt) - 












9.8 
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Figure 9-l» Analog Computer Program for Fourth Order Fade Approxiinant, 




1^1 



TABLE OF COMPUTER COMPONENTS 



Potentiometer Settings 


Component Values 


®i “ '^i^i 


1 Megohm 


32 = R2CJ, 




33 = R3C2 


'3 s 


= 0.2 C2 


o^ 

00 




R9 Rii 


^ = V3 


^12 ^14 


a^ = 0,467 C3 


^f6 



*9 ■ V4 


2 Megohm 


"10 ■ ‘=4 


^f5 


■ ''u ‘^4 
R* « 


1 Microfarad 






R,/ 

"f6 


"3 "4 



Th<? time is scaled such that computer time is equal 
to ten times problem time. 
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